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Radar  Cross-Section  of  Damped  Cylinders 
and  Anisotropic  Cylinders 


This  quarter,  we  have  generalized  the  TDFD  RCS  code  to  handle  materials 
characterized  by  ■, 


(''1)  Lossy,  frequency -dependent  dielectrics  } 
Cl)  Magnetic  loss  and  conductivity  j 
(3)  Anisotropy  in  the  xy  plane. 

(U)  Embedded  resistive  or  reactive  sheets  , 


Our  treatment  of  frequency  dependence  goes  beyond  what  had  been  the 
state  of  the  art  in  TDFD  codes;  the  frequency  dependence  here  may  be 
anisotropic  and  may  extend  across  a  molecular  resonance.  The  techniques  for 
doing  this  are  documented  in  the  attached  article,  "Time-domain  treatment  of 
Maxwell's  equations  in  frequency  dependent  media.*  This  article  was  written 
for  the  March  conference  in  Monterey  and  represents  work  entirely  inspired 

by  the  present  contract.  /, 
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A  second  computer  code  has  also  been  written  to  compute  the  RCS  of  a 
circular  cylinder  composed  of  an  arbitrary  number  of  concentric  shells,  each 

ic 

characterized  by  an  arbitrary  a,  a  ,  e  and  \i.  This  code  uses  the  expansion 
of  plane  wave  in  cylindrical  harmonics  and  matching  of  harmonic  coefficients 
at  each  interface.  It  is  described  in  the  attached  note,  "Scattering  from  a 
layered  dielectric  cylinder." 


Code-code  comparisons  between  the  TDFD  and  cylindrical  harmonic  codes 

have  been  run  for  a  perfect  conductor  .5  m  in  radius,  bare  and  covered  by  a 

★ 

damper  .5  m  thick.  The  damper  has  properties  «/«„  -  -  1;  a  -  eo  /n  -  4 

x  10 ’■*.  These  values  were  selected  to  give  a  skin  depth  on  the  order  of  the 
damper  thickness  at  frequencies  (50  -  500  MHz)  for  which  calculations  were 
run.  The  TDFD  code  utilized  square  cells  4  cm  on  each  side  or  25  cells  to  a 
cylinder  diameter. 
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While  our  TDFD  code  only  treats  the  TM  case,  we  are  able  to  simulate 
the  TE  problem  through  appeal  to  duality.  In  particular,  the  previously 
described  code -code  comparison  was  rerun  with  the  cylinder  perfectly  mag- 

iMr 

netically  conducting  (o  —  ®,  a  —  0) ;  and  the  damper  unchanged. 

The  following  four  figures  present  code-code  results  for  the  four 
possible  combinations  of  electric  and  magnetic  conducting  cylinders  with  and 
without  damping  shells . 
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Scattering  from  a  Layered  Dielectric  Cylinder 


This  note  discusses  scattering  of  a  plane  wave  by  a  circular  dielectric 
cylinder  composed  of  concentric  layers  of  different  materials.  Let  us 

lit 

assume  there  are  N  layers,  with  layer  i  characterized  by  e^,  (i a a .  and 
outer  radius  a^ .  We  shall  here  treat  the  TM  case  (H  along  the  cylinder 
axis;  E  transverse),  although  the  TM  problem  is  nearly  identical 
mathematically . 

Assume  the  incident  wave  is  propagating  in  the  +  y  direction. 


HinC(£,t)  -  Y0ize 


i(k0rsin^  -  wt) 


.  i(k„rsin*  -  ut) 

ElnC(£.t)  -  -  j^e 


Here,  Y0  t^ie  admittance  of  free  space.  The  same  symbol  will  sub¬ 
sequently  be  used  to  designate  Neuman  functions,  but  context  should  keep  the 
meaning  unambiguous.  Additionally,  k0  is  the  free  space  wavenumber,  uje 0n0 
—  w/c . 


Equation  (1)  may  be  expanded  in  cylindrical  harmonics, 


Hinc(r,t)  -  Y0Iz  ^  Jn(k0r)ein* 


-  Y0i  [j0(k0r)  +  'S  2J  (k0r)cos  n<£  +  S  2iJ  (k0r)sin  n<$] 


In  the  future,  it  will  be  useful  to  designate  the  coefficients  of  these 

harmonics  as  a^nc; 
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aLnc  -  2i 
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n  >  0,  even 


n  odd 


Since  V  x  H 
becomes 
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jwf0E  ,  the  cylindrical  harmonic  expansion  for  E 


iY  i  40  “ 
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The  innermost  material  will  include  the  cylinder  axis.  Thus,  only 
Bessel  functions  of  the  first  kind  are  permitted  in  the  solution  there: 


H1 (it  t) 


00  00  -| 

Yliz  ^  aiJn(k1r)cos  n*  +  ^  anJn(k  ir)sin  n<^ 
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a'J' (k,r)sin 
n  n  1 


(7) 


where  Yx  is  the  admittance  of  medium  1 


+  jt^/w 
•k 

+  j a^/u 


(8) 


and  kx  is  the  wavenumber  of  medium  1, 


ki 


+  j^i/w)(^i  +  j^i/w 


(9) 


The  N-l  concentric  shells  will  permit  solutions  of  both  kinds.  Thus, 
in  region  i,  1  <  i  £  N,  a^J  (ktr)  of  eqs.  (6) -(9)  becomes  replaced  by 


a1 J  (k,r)  - 
n  n  1 


a1!  (k.r)  +  bS  (k-r) 
n  nv  i  '  n  n  i 


(10) 


Finally,  in  free  space  outside  the  cylinder,  the  first  Hankel  function 
is  the  only  permitted  solution  for  the  scattered  field.  Thus,  in  this 

inC 

region,  a^  Jn(k0r)  is  replaced  by 
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scat..(l).,  . 
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n 
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(11) 


in  eqs .  (3) - (5) . 


For  each  azimuchal  harmonic,  eqs.  (12) -(14)  comprise  a  set  of  2N  linear 
equations  in  2N  unknowns.  The  associated  matrix  is  five-banded,  and  ex¬ 
tremely  easy  to  solve  by  Gaussian  elimination.  (The  main  diagonal  and  first 
diagonal  off  each  side  of  the  main  is  full.  The  second  diagonal  off  each 
side  of  the  main  is  half  zeros.) 

The  quantities  of  interest  in  RCS  evaluation  are  the  ascat's.  A  two- 

n 

dimensional  bislatic  RCS  is  defined  by 


RCS(^)  -  2w  r^1"® 


Einc 


(15) 


The  scattered  electric  field  (neglecting  the  reactive  radial  component)  is 


Escat 


(£.t) 


t  °° 


^  ascatH^l)  (fcor)cos 
n-0,2 


a®catH;(k0r)sin  n 4 

n-1,2 


(16) 


Here,  use  is  made  of  the  identity  k0Y0  -  we0. 

For  large  arguments ,  Hankel  functions  have  the  asymptotic  limit 


(17) 


Substitution  of  eqs.  (16)  and  (17)  in  (15)  then  yields  the  RCS  in  terms  of 
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ABSTRACT 

Longmire  and  Longley  (DNA  report  3167F,  1973)  have  described  a  method 
for  representing  certain  types  of  frequency -dependent  media  in  time-domain 
finite  difference  codes .  They  assumed  the  media  could  be  described  by  a 
series  of  relaxation  phenomena,  one  phenomenon  per  decade.  Each  phenomenon 
contributes  a  "current"  to  the  7  x  H  equation,  although  these  "currents"  do 
not  conveniently  fit  either  the  concept  of  conduction  or  displacement. 

Recent  discussions  with  ONR  personnel  have  raised  interest  in  more 
general  types  of  media  characterization.  Specifically,  there  is  now  a 
desire  to  input  the  SEM  or  Prony  parameters  of  a  medium’s  o  and  e  to  time- 
domain  codes.  In  other  words,  we  now  wish  to  work  with  the  actual  poles  and 
residues  of  a  and  e  rather  than  assuming  the  poles  are  uniformly  spaced  one 
per  decade.  The  present  work  describes  how  this  may  be  done. 

It  also  goes  beyond  the  work  of  Longmire  and  Longley  in  permitting  use 
of  complex  as  well  as  real  poles.  Prony  "currents"  associated  with  complex 
poles  obey  a  temporally  second-order  differential  equation,  as  opposed  to 
the  Prony  "currents"  of  real  poles.  In  both  this  formulation,  and  of  the 
Longmire -Longley  formulation,  the  latter  obey  temporally  first-order  dif¬ 
ferential  equations. 

In  actual  time-domain  finite-difference  codes,  it  is  most  exact  to 
evaluate  the  Prony  "currents"  as  well  as  the  electric  fields  at  the  same 
spatial  and  temporal  points.  In  the  past,  this  simultaneous  solution  has 
not  been  implemented,  due  to  the  peculiar  coupling  of  the  equations  for 
electric  fields  and  Prony  "currents".  For  the  case  of  all  poles  real,  we 
here  also  describe  how  to  perform  this  exact  solution  using  state  theory; 

Work  was  sponsored  by  Office  of  Naval  Research  and  DARPA  under  Contract  No. 
N00014-85-C-0852 
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i.e.,  by  exponentiating  the  matrices  coupling  all  the  electric  fields  and 
Prony  "currents". 


INTRODUCTION 

Consider  a  medium  with  anisotropic,  frequency -dependent  electrical 
properties.  The  electrical  response  of  such  a  material  may  be  fairly 
generally  described  by 


V  x  H  -  J(t)  +  Jf(t) 


(1) 


ii 


where  J^t)  is  a  forced  current  and 


V 


« 


!(t)  -  a0-  E(t)  +  t"  '  D£(t)  +  K(t-t')  ‘  D£(t')dt'  (2) 


with  D  indicating  the  time -derivative  operator.  In  this  formulation,  o0, 
«  ,  and  K(t-t')  are  second- rank  tensors. 

MO  ■ 

The  frequency -domain  form  of  eq.  (2)  is 


t; 


ji(«)  -  £a0  +  iwf  r  +  iw  J^e‘iwuK(u)du]  •  £(«) 


(3) 


i 


Separation  of  eq.  (3)  into  real  and  imaginary  parts  gives  representations 
for  the  frequency -dependent  conductivity  and  permittivity  tensors, 


o<w)  -  £0  + 


Re  iwj^e " iwuK(u) duj 


(4) 


to 

V-. 


14 


«<«)  -  +  Re^j7e'iwUK(u)duj 


(5) 


Longmire  and  Longley^  have  considered  the  scalar  version  of  this  for¬ 
mulation  for  the  special  case  when  K(u)  can  be  expressed  as  an  exponential 
series , 


M 

\  ~P  u 

K(u)  -  /  a  e  m  (6) 

m-1 


For  this  situation,  eq.  (2)  may  be  rewritten 


I(t) 


ffoKt)  +  «.!>£(  t)  + 


M 

Lv.- 


(7) 


with 


vc) 


(8) 


Equation  (8)  is  equivalent  to  the  differential  equation 


°4<t)  -  D£<c>  - 


(9) 


Longmire  and  Longley  assumed  that  materials  could  be  represented  by  the 
exponential  series  of  eq.  (6)  with  one  term  for  each  decade  of  frequency 
over  the  spectrum  of  interest.  This  is  equivalent  to  doing  a  Prony  expan¬ 
sion  of  K(u)  (or  ff(u)  and  «(«)]  with  the  poles  forced  to  be  spaced  at 


While  this  assumption  has  been  claimed  to  be  reasonably  accurate  for  wet 
soil,  it  would  seem  generally  more  correct  to  determine  the  poles  from  a 
Prony  analysis  of  the  medium's  measured  frequency -dependent  characteristics. 
This  is  especially  likely  to  be  true  if  the  material  exhibits  rapid  varia¬ 
tion  in  a  and  e  with  frequency. 

STATE  THEORY  APPLICATIONS 

Let  us  first  asstime  the  Prony  analysis  reveals  no  complex -conjugate 
pole  pairs.  In  general,  the  am  will  be  second  rank  tensors,  but  the  will 
only  be  scalars.  Then  for  every  pole,  each  component  of  J  will  obey 


DEi  *  wmi  *  Vmi  “  0  »  "  1  *  M.  i  -  1  -  3 


(ID 


Additionally,  the  tensor  form  of  eq.  (7)  gives 


M 


.  .DE. 

j 


+  o. 


•i1E1  +  ^  J 
J  m-1 


mij  mj 


J  ,  -  J  j 


i,J  -  1  *  3 


(12) 


where 


Jj  -  (V  x  H  -  if)j 


(13) 


Equations  (11)  and  (12)  constitute  a  set  of  3(M+1)  coupled  first  order 
differential  equations. 


If  anisotropy  and  frequency  dependence  were  not  present,  the  usual 

method  of  numerical  solution  would  be  explicit  time-domain  finite 

differencing.  In  this  method,  £  and  H  evaluation  points  alternate  both 
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spatially  and  temporally  using  a  well- tested  leapfrog  arrangement.  In 
this  arrangement,  no  two  equations  are  coupled,  and  e”+^^(I, J,K)  means  Ex 
evaluated  at  ((I  +  1/2)AX,  JAY,  KAZ,  (n  +  l/2)At). 

However,  the  present  system  of  equations  requires  the  three  Ej's  and  3M 
J^'s  all  to  be  evaluated  simultaneously.  While  this  cannot  be  done  using 
conventional  time-domain  finite  differencing,  state  theory  does  indicate  an 
appropriate  generalization  of  time-domain  finite  differencing. 

First,  let  us  consider  the  case  where  anisotropy,  but  not  frequency 
dependence ,  is  present , 


[«J  D[E]  +  [o0][ E]  -  [J] 

This  matrix  differential  equation  has  a  homogeneous  solution 


(14) 


[E]h 


-Uj'VoJt 


[A] 


(15) 


and  a  particular  solution 


[E]p  -  [tf0] ~ 1 t J ] 


(16) 


giving  a  general  solution 


IE]  - 


K 


[A]  +  [*0]'*[J] 


(17) 


The  constant  vector  [A]  may  be  evaluated  at  (n  -  1/2) At: 


[E]n'1/2  -  [A]  +  [*0]-1[J]n 


(18) 


This  gives  the  new  E-field  vector  in  terms  of  the  old, 


[E] 


n+1/2  _ 


-[«, 


yh*  oiAt 


[E] 


n-1/2 


♦  t- 


-[«.]’ Vol  At! 


J[* 


0]'1[*nTl 


(19) 


Similar  exponential  matrix  techniques  have  been  reported  for  time- 
domain  solution  of  generalized  multi-conductor  transmission  lines. ^  In  the 
previous  work,  one  may  see  how  to  evaluate  eq.  (19)  if  [o0]  is  singular  or 
if  [ e^] *^[a0]At  has  arbitrarily  large  elements.  Basically,  matrices  are 
exponentiated  using  the  power-series  representation  of  an  exponential. 

If  frequency  dependence  is  present,  the  [E]  vector  of  eqs.  (14) -(19) 
becomes  replaced  by 


[E]  - 


-  [E'J 


The  [e  1  matrix  becomes 

flO 


(20) 


(21) 


and  the  [o0]  matrix  becomes 


2m 


[<*<>] 


°0 


2i 


•Pi i  *•*  o 


-  [*'] 


(22) 


Lastly,  the  forcing  vector  becomes 


[J] 


1 

Q. 

fi 


U'l 


(23) 


Then  the  matrix  equation  for  simultaneous  advancement  of  £  and  the  is 
[E']n+1/2  -  e'1*^  [ff']AtjE, jH*1/2  +  .  e-t€'l  [ff']dtJ[ff,j.ljj,jn  (24) 


In  the  past,  time-domain  finite  differencing  has  not  often  considered 
anisotropy.  Frequency -dependent  effects  have  been  included  by  using  the  old 
J^"^2  to  find  the  new  fi1**^2.  (This  decouples  £  from  the  ^  in  eq.  (12).) 
Then  the  new  £n+^2  have  been  used  to  find  the  new  J^+^2  from  eq.  (11). 


TREATMENT  OF  COMPLEX  POLES 

If  Prony  analysis  of  the  material's  frequency  dependence  reveals  com¬ 
plex  pole  pairs,  a  more  general  treatment  becomes  necessary.  In  this  case, 
K(u)  will  contain  terms  of  the  form 


b  sin(y  u  +  d  )e 
— n  ’m  m 


-  b  cosd  sin7  ue  m  +  b  sini  C0S7  ue 
— m  m  m  -m  m  m 


The  J  (t)  of  eq.  (8)  now  becomes 


pt  -p 

J  (t)  -  J  D£(t')(cosi  sin7  (t-t')e 


+  sin^  cosy  (t-t')e 

m  m 


“  ^c(t)cos*o  +  4s(t)sin*m 


where  and  are  the  parts  of  associated  with  cos^m  and  s in^m , 

respectively: 


ft  -P  (t-t') 

4c(t)  ”  LD£<t')sin7m(t-t')e 


ft  -Pm<t-t'> 

4s(t,)  "  J.Dfi(t')cosVt't')e  dt 


Differentiation  of  J  (t)  and  J^ft)  yields 


D^c<C>  “  -^c<C>  +  Vms<C> 


IS? 
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wytvvvjiraiw  ,wjww  *g  lpjvj^j  wj  r/Tjrj™/w  ■'  ■  ■ 


W»s<'>  -  DECt)  -  V„(t)  1-  7rfimc<t> 


(30) 


These  equations  can  be  solved  for  J  (t)  and  J|||g (t)  using  Heavyside  algebra: 


[<»  +  />„>’  *  7j]  4c<t>  -  7„DE<t)  (31) 

[(D  +  0ay  7  7*1  J^ft)  -  (D  +  ^)DE(t)  (32) 


Thus,  J  (t)  obeys  the  differential  equation 


[(D  +  fia)*  +  7*1  ^(t)  -  [sin^m(D  +  flj  +  cos^lDKt)  (33) 


In  principle,  equations  like  (33)  could  be  added  to  the  set  of  equa¬ 
tions  given  by  (11)  and  (12),  and  the  entire  ensemble  solved  by  state 
theory.  This  approach,  however,  requires  treatment  of  second-order  matrix 
differential  equations  of  the  form 


[A]D*[E]  +  [B)D[E]  +  [C}[E]  -  [F] 


(34) 


The  homogeneous  solution  of  this  equations  includes  square  roots  and  complex 
exponents  of  matrices;  it  is  much  more  complicated  than  eqs.  (14)  -  (19). 
(To  the  best  of  our  knowledge,  exponential  differencing  has  never  been 
applied  even  to  scalar  second-order  differential  equations.) 

Consequently,  when  Prony  analysis  of  the  material  data  yields  complex 
pole  pairs ,  our  present  strategy  is  to  fall  back  to  the  old  technique  for 
dealing  with  real  poles:  First  find  using  the  old  .  Then  use 


REFERENCES 


Longmire,  C.L.  and  H.J.  Longley,  "Time  Domain  Treatment  of  Media  with 
Frequency -Dependent  Electrical  Parameters,"  DNA  3167F,  23  September 
1973. 

Yee,  K.S.,  "Numerical  Solution  of  Initial  Boundary  Value  Problems 
Involving  Maxwell's  Equations  in  Isotropic  Media,"  IEEE  Trans.  Antennas 
&  Pronaeat. ,  Vol.  AP-14,  pp.  302-307,  May  1966 

Longmire,  C.L. ,  "State  of  the  Art  in  IEMP  and  SGEMP  Calculations,"  IEEE 
Trans.  Nuc.  Sci. .  Vol.  NS -22,  pp.  2340-2344,  December  1975. 

Holland,  R. ,  "THREDE:  A  Free-Field  EMP  Coupling  and  Scattering  Code," 
IEEE  Trans.  Nuc.  Sci. .  Vol.  NS-24,  pp.  2416-2421,  December  1977. 

Holland,  R. ,  "The  EMP  Tunnel  Code  DHTRANS,"  IEEE  Trans.  Nuc.  Sci..  Vol. 
NS -25,  pp.  1404-1411,  December  1978. 


